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Bloch equations for the atomic population and the polarization/coherence and the equation of 
motion for the photon number in a laser are solved in steady state as a function of the pump rate. 
Two level atom and two modes of three levels atom are investigated. Close to threshold the usual 
linear dependence of the intensity on the pump rate in found for all cases. However, far above 
threshold strongly nonlinear dependence is encountered. In the cases for which the pump connects 
the lower lasing state to one of the excited states the character of the non-linearity differs crucially 
from the cases when the pump in not related directly to the lower lasing state. Non-monotonic 
dependence of laser intensity upon the pump rate is predicted. Detailed discussion of the nonlinear 
behavior is presented, including saturation and depletion effects. 
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The behavior of the number of photons and the laser intensity as function of the pump rate near and above 
threshold has recently been studied [0. In particular, linear dependence of the intensity close to threshold and some 
non-linearity above threshold has been found [J. The purpose of this paper is to point out the non-linearity of 
the laser intensity as function of pump and the crucial role of the rate of depletion of the lower lasing state in the 
dynamics of lasers. It is usually believed [Q that the faster is the rate of depletion of the lower laser level, the better, 
because this gives rise to a larger inversion at the lasing transition. However, as it is shown in this paper, there is a 
limitation on the depletion rate of the lower lasing level. For any arbitrarily large pump rate there exists a critical 
' value of the depletion rate of the lower lasing level, at which the lasing is broken down. 

in 

In this paper we study closed two and three level systems with incoherent pump. Consider first two level scheme 
shown on Fig. 1. The corresponding optical Bloch equations are 
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Au = -7/>ii + r /°oo - ig(zpoi - c.c), (1) 
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pio = -TiPio + igz(pu - poo), (2) 
h = —2nn + iNg(zpoi — c.c), (3) 

where N is the total number of atoms, g is the electric dipole coupling, k is cavity decay rate, n — z*z is the number 
of photons inside the cavity, 7 and T are spontaneous emission rate and the pump rate respectively. The steady state 
d ' solution for the photon number n reads 
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where P = T/j is the relative pump rate. It follows from eq. (Q) that in order to get lasing the following condition 
must be satisfied 
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P>P thr = ^£, (5) 
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which is usually referred to as a threshold condition. The necessary condition for eq. (0) to be true is 

7, < ^. (6) 

K 

If this inequality is not satisfied, lasing cannot be obtained at any pump rate. In case of Ng 2 / 'ivy± 3> 1, when eq. 
(||) is satisfied automatically, the threshold condition is simple 

P > 1 (7) 
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which means that the pump rate T should exceed the decay rate 7 of the upper lasing level. Inequality gives 
rise to a restriction on the pump rate. Indeed, it follows from the derivation of the relaxation rates (see Refs. 0)) 
that 

r + 7 + i P h m 

where "f p h stands for the contribution from the phase destruction processes like elastic collisions. Combining eqs. 
@" @ we obtain the following restriction for the pump rate P 

1 <P< 2 J^-1-S, 5=3*. (9) 
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At hrst glance at eq. (^) , the dependence of the mean photon number n upon the pump parameter P seems to be 
linear. However this is not true due to the dependence of the transversal relaxation rate 7^ upon the pump parameter. 
Indeed, taking into account eq. (0) we get, instead of eq. (0), the following expression for the mean photon number 



^-\P - 1 - (P + 1)(P + 1 + 5)-^ 
4k 1 >2Ng 



P-l-(P+l)(P+l + <5)^ 7 ^]. (10) 



The dependence (|10|) is not linear, it is rather quadratic in P (see Fig. 2 ) with maximum at 

Po = -^--l--. 11 

kj 2 

Note that condition ([}]) is necessary but not enough to get lasing. The region of the pump parameter P in which 
lasing is really possible is given by (see fig. 2) 

KP<^-*-3, (12) 

Here we assumed that Ng 2 /kj 3> 1 which is normally true. Maximum value of the photon number is defined by 



eqs. (p. 0D, flliP and has the following form, provided Ng 2 /n~/ » 1 

n max « l{^) 2 (13) 



The non-linearity in eq. (10) stems from the fact that according to eq. (||) the rate of transversal relaxation j± 
depends upon the pump rate T. The larger is the pump rate, the shorter is atomic radiative life time which is inverse 
proportional to "fperp. Thus, a "too fast" pump can give rise to a reduction of the photon number (or intensity) and 
can even switch off the lasing process. At relatively low pump this reduction is compensated by the gain which is 
given by the linear term in eq. (|l0|). However, when the pump parameter P becomes large enough, i.e. of order of 
P ~ Ng 2 /n~/, the non-linearity plays an important role and cannot be neglected. The two level model exhibits an 
unusual non-linearity with respect to pump rate. 



To get a more realistic model one has to involve multilevel schemes. In such schemes however, the depletion of the 
lower lasing level is not necessarily associated with pumping mechanism and therefore, the character of non-linearity 
may be quite different. 

Let us consider now a more realistic three level model shown in fig. 3, which is a particular case of a general three 
level model considered in Ref. ||. Note that the two schemes depicted in figs. 3a and 3b are equivalent mathemat- 
ically, i.e., they are formally described by the same set of equations. However, physically the dynamic properties of 
these two schemes are completely different because of the different ways of pumping as described by the role of the 
relevant parameters (see figs. 3a and 3b). In the scheme of fig. 3a the lower lasing state |0 > is depleted with the rate 
702 to the ground state 2 >, and then the pumping is used to excite the atom to the upper lasing state |1 > with 
pump rate 721. In the scheme of fig. 3b the pump mechanism is used to excite the atom from the lower lasing state 
> which is now the ground state, to the upper lasing state |2 > which is then depleted with rate 721. In the last 
case the role of the pump rate is played by the rate 702 . As will be seen in the following this difference in the pump 
mechanism is crucial for the dynamics of laser. 
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Let us write down a set of optical Bloch equations for the three level model (fig. 3). 



pn = 721P22 - 710P11 - g(izpoi + c.c), (14) 

(>22 = 702P00 - 721P22, (15) 

pio = -7-LPio + igz(pn - Poo), (16) 

h = —2nn + Ng(izpoi + c.c.) (17) 



The steady state solution for the number of photons inside the cavity reads 

N 721(702 — 710) 7-1 702721 + 702710 + 721710 

n = r . 

2k 702 + 2721 2g z 702 + 2721 



(18) 



This formulae is valid for both models shown in figs. 3a and 3b. It seems from eq. (18) that the photon number n 
depends upon both rates 702 and 721 in a similar way, i.e. both dependencies are of linear-fractional type, which means 
that n tends asymptotically to some fix value when increasing either the rate 702 or 721. This type of non-linearity 
has been brought about by the nonlinear dependence of both longitudinal relaxation rate and equilibrium inversion 
A upon the rates 702 and 721. This can be shown explicitly by rewriting eq. (|q) in the following form 
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where 



^7|| A 7X7|| (ig) 



_ 2(721702 + 702710 + 721721) ,„„>. 

7|| = , (20) 

702 + 2721 

A _ 721(702 - 7io) / 21 s 



721702 + 702710 + 721721 



However, as it has already been mentioned in the discussion of the two level model, the dipole relaxation rate 7j_ 
is related to the other rates. In the three level case under consideration the following relation for j±_ takes place 

7i = ^(Tio + 702 + Iph), (22) 

where j p h stands for collisional dephasing rate. Note that j± depends on 702 but does not depend upon the rate 
721. It is this asymmetry that is responsible for the difference in dynamic properties of lasers with the above two 
different pump mechanisms. Indeed, in the case of the scheme of fig. 3a, 7j_ does not depend upon the pump rate 721 
and hence the dependence of the photon number upon pump remains linear-fractional. In contrast, in the case of fig. 
3b the role of pump rate plays 702 and thus 7j_ brings about additional dependence upon the pump rate. Therefore 
the type of non-linearity of the photon number as function of pump in this case differs from that of the scheme 3a. 
This difference is crucial for the dynamics of the lasers with the different pump schemes. 

Let us consider the two above mentioned modes of pumping separately. In the case of the scheme on fig. 3a 



equations (18) and (E2J) result in the following expression for the photon number 



where 



A Pi(l - ei) - + ei + fr)[Pi(l + 61) + ei] 
11 l + 2Pi [ ' 

_ N^/ 2 _ K702 _ 710 D _ 721 , _ l P h foA , 

Ai = — — , si = ——5, ei = , Pi = , St = (24) 

2k 2Ng z 702 702 702 

In fig. 4a rii is plotted as function of pump parameter P±. At the beginning, when Pi << 1, n\ grows linearly 
with the pump rate Pi. However, when the pump is strong enough, i.e. 721 > ^S 2 -, the dependence ni(Pi) becomes 
non-linear. The origin of this non-linearity lies in the large disparity between the rates of depletion of the states |2 > 
and |0 >. When the lower lasing state |2 > is depleted much faster than the ground state |0 >, there is a bottleneck at 
the transition |0 > —\2 > so that further increase of the pump rate does not result in increase of the photon number. 
This kind of non-linearity was mentioned in Ref. || . 
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Threshold condition for the scheme 3a is 



P > ei.si(l + ei + £i) 

1 l-ei-.si(l + ei+<5i)(l + ei) 1 ' 

It follows from eq. (|2^) that there are two additional necessary conditions for getting lasing: 

ei < 1, (26) 

51 < + + (2?) 

As it follows from the definition of the equilibrium inversion, A, (see eq.(pl|) ), the condition fl26| ) merely means 
that there should be initial positive inversion at the lasing transition. As for cq.(p7|) we can rewrite it assumming 
that 2Ng 2 << 1 m the following form 

1<^<^!_^_3 (28) 
7io K710 710 

which is similar to eq. (^2[) . We came to the same restriction for the rate of depletion of the lower lasing level as 
for the two level model. This time however, this does not result in any restriction for the pump rate 721, which can 
be readily understood because in contrast to the two level scheme, in the scheme under consideration, there is no 
connection between the transversal relaxation rate and the pump rate. Inequality (^) has also another important 
meaning which becomes obvious if we rewrite it in the following way 

^V>iV^ = ^ (l + ei 1 +W + £l) . (29) 
2g 2 1 - ei 

Thus there is a minimal number of atoms N m i n needed to get lasing. If N < N m i n one cannot get lasing at any 
pump rate. 



Let us consider the scheme shown in fig. 3b. Now instead of eq.(|23j) we have for the photon number 

. P-2 - e 2 - s 2 {P 2 + e 2 + S 2 ){P 2 + e 2 + P 2 e 2 ) , on , 

«2 = A 2 , (30) 

where 

\ _ ^721 _ K721 _ 710 _ 702 x _l P h , Q1 , 

A 2 = -z — . s 2 = — T , e 2 = , P 2 = , 6 2 = — ■ (31) 

2k 2Ng 2 721 721 721 

Note that the third term in the numerator contains additional dependence upon the relative pump rate P 2 . This 
additional dependence has been brought about by the fact that the dipole relaxation rate 7^ depends on the pump 
rate 702- The photon number n 2 is plotted in fig. 4b as a function of the pump rate. Again, as in the previous case, 
the photon number grows linearly at relatively weak pump and becomes non-linear at strong pump. However, the 
character of this non-linearity is crucially different. It has the same form as in the case of the two level scheme (cf. 
fig. 2). This difference is due to the above mentioned additional dependence of the photon number upon the pump. 
The photon number has a maximum at 

Ng 2 S 2 , x 

P2 = — - i - e 2 . 32 
K721 2 

Here we assummed that e 2 << 1. Eq.(|32]) is similar to eq.([ll|) for the two level model. In order to get lasing the 
following double inequality must be satisfied, provided S2> e 2 << 1 

1< P 2 < — - S 2 (33) 
s 2 

which is similar to eq.(|l^) for the two level model. Again we got the same restriction for the rate of depletion of 



the lower lasing state |0 > as in the case of eq.(28). This time however, this is the restriction for the pump rate P 2 



since when P 2 reaches the value of l/s 2 — S 2 , the lasing is broken down for the same reason as in the case of two level 
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model. 



In summary, we have shown that the nonlinear dependence of the photon number upon the pump parameter is very 
sensitive to the type of pump mechanism. When the pumping is used to excite the lower lasing state to one of the 
upper atomic states, the laser exhibits a peculiar behavior shown in figs. 2 and 4b for two and three level schemes 
respectively. The number of photons as a function of pump reaches its maximum and then slows down untill some 
critical point at which the lasing ceases. It should be noted that the described non-linearity manifests itself at very 
large pump rates when the ratio p/pthr reaches a few orders of magnitude. Such a regime may be difficult to achieve 
in conventional lasers. For instance, in order to reach the non-linear region in experiments with GaAs laser used in 
Ref. the ratio p/pthr should be of order of 10 5 . Nevertheless, it may be realized in a microlaser. 
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Figure captions 



1. Two level scheme. 

2. The photon number vs relative pump rate for the two level scheme. The parameters are A = 10 3 , S = 10 5 , s = 
7 1(T 7 , 1(T 6 and 1.33 1CT 6 for curves 1, 2 and 3, respectively. 

3. Three level schemes with different pump mechanism. 

4. Photon number vs pump rate for the pump schemes shown in figs. 3a and 3b. The parameters are: Ai = 10 6 , 
ei = 0.01, <5i = 0, si = (curve 1), 0.2 (curve 2) and 0.5 (curve 3) for fig 4a; A 2 = 10 5 , e 2 = 0, S 2 = 0.1, s 2 = 0.1 
(curve 1), 0.02 (curve 2) and 0.01 (curve 3) for fig 4b. 
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